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Abstract - We discuss the 2+1 flavor Polyakov loop enhanced Nambu-Jona-Lasinio model in a 
finite volume. The main objective is to check the volume scaling of thermodynamic observables 
for various temperatures and chemical potentials. We observe the possible violation of the scaling 
with system size in a considerable window along the whole transition region in the T — p q plane. 
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-d! Introduction. — The hot and/or dense matter cre¬ 
ated in ultra-relativistic heavy ion collisions is supposed to 
(N possess a rich phase structure. In the intermediate regime 
of temperature and baryon chemical potential in the range 
in of few hundred MeV, the defining characters of the phases 
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are the color confinement and chiral properties. While for 


qq low baryon densities the matter has a smooth crossover 
c from color confined chiral symmetry broken phase to color 
deconfined chiral symmetry restored phase, at high enough 
densities, this transition may be of first order. A critical 
l/~) end point seems to naturally occur in such a situation. 
Establishing this scenario forms an integral part of explo- 
ration in the international collaborative experiments at 
•’“j CERN and BNL and the upcoming experiments at GSI. 
rj| The matter formed in heavy-ion collision experiments 
has a finite volume, which depends on the size of the col¬ 
liding nuclei as well as the center of mass energy (y/s) and 
the centrality of collisions. Measurement of this system 
size is quite non-trivial. There are several estimates of 
the system size at freeze-out for different y/s and different 
centralities from the measurement of HBT radii |T], which 
indicate that the freeze out volume increases as the y/s 
increases. The freeze out volume was found to be in the 
range of 2000 fm 3 to 3000 fm 3 . On the other hand in the 
URQMD model [2] the volume of homogeneity has been 
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calculated and compared with the experimentally avail¬ 
able results in Ref. [3]- The homogeneity volume mu is 
given by the product of the HBT radii R ou t, Rside and 
Riong , signifying the fraction of the fireball from which 
particles in a particular momentum window is emitted. 
The volume of homogeneity was obtained in the transverse 
momentum window of 300-400 MeV at different centrali¬ 
ties and varying y/s from 62.4 GeV to 2760 GeV. These 
UrQMD simulations show that the homogeneity volume 
has a power law scaling with multiplicity and varies in the 
range of 50 to 250 fm 3 . The ALICE collaboration has 
also estimated the freeze-out volume to be in this range 
for different y/s and different colliding nuclei [6]. Given 
that these are the freeze-out volumes, the initial volumes 
are expected to be much smaller m- 

Theoretical developments to understand the effects of 
finite volume on the strongly interacting matter is repre¬ 
sented in Refs. EH2D- General overview of the develop¬ 
ments in this direction is reviewed in some of our recent 
work [32lf34l . In case of high energy heavy-ion collisions it 
was shown in Ref. [35] that the finite size effects give rise 
to important consequences in the later stage of evolution 
of the hadronic bubbles. In Ref. [55], the authors have 
discussed the importance of finite size scaling in order to 
properly locate the critical end point (CEP) on the phase 
diagram and emphasized on its use for the experimental 
data analysis program EZj- In Ref. [32] some of us have 
studied the thermodynamic properties of strongly inter¬ 
acting matter in a finite volume using Polyakov-Nambu- 
Jona- Lasinio (PNJL) model for 2 and 2+1 flavors. It was 
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found that the cross-over temperature at zero baryon den¬ 
sity decreases with decrease in volume. Furthermore the 
critical end point at finite temperature T and chemical 
potential p goes towards higher p and lower T domain as 
system size is reduced. For a system confined to a lateral 
size R = 2/m, it was found that the CEP vanishes and the 
whole phase boundary becomes a cross-over line. We also 
discussed the possibility of chiral symmetry restoration in 
a color confined state for a system with finite size. 

Fluctuations of conserved charges are sensitive indica¬ 
tors of the transition from hadronic matter to partonic 
state. Also the existence of the CEP may be signalled 
by the diverging behavior of fluctuations. From the chiral 
susceptibilities in the 2-flavor quark-meson model using a 
renormalisation group approach [381 the CEP was found 
to have lower temperature and higher chemical poten¬ 
tials with decreasing volume, similar to our results in the 
PNJL model [32] • Similar agreement was reported from 
baryon number susceptibility ratios via Dyson-Schwinger 
approach using 2-flavor system [3[£. On the lattice, for a 
pure gluon theory, the Polyakov loop susceptibility has 
been calculated for a finite volume [32 showing that 
the transition temperature may increase with decreasing 
system size. In a lattice version of the 2-flavor PNJL 
model [40] the diagonal quark number susceptibility was 
found to decrease with increasing volume. Our recent 
studies in the mean-field approximation of 2-flavor PNJL 
model [33] as well as in hadron resonance gas model [341 . 
give similar results. Additionally we found that the vol¬ 
ume scaling of susceptibility ratios is violated near the 
cross-over region in the PNJL model and near the freeze- 
out surface in the hadronic model. It is therefore im¬ 
portant to explore the situation in the 2+1 flavor PNJL 
model, which is reported in the present work. 

We first give a brief description of the PNJL model fol¬ 
lowed by the results for the fluctuations and correlations. 

Model. — The NJL model BTffllJI gives a satisfactory 
description of strongly interacting matter at zero tempera¬ 
ture and chemical potentials. However this is a model with 
global color conservation and therefore no confinement of 
color charges appear at non-zero temperatures and den¬ 
sities. A suitable modification is done by introducing a 
background field that mimics the behavior of the Polyakov 
loop to bring in the effects of confinement BZH55]. Consid¬ 
erable progress has been made to have an understanding 
of strongly interacting matter using the PNJL model (see 
e.g. 32 3311+11 [731 . Here we have considered 2+1 flavor 
PNJL model including upto six-quark interactions Ena. 
To include the effect of finite volumes, a non-zero lower 
momentum cutoff p m in = 7r/i? = A is introduced. Here 
R is the lateral size of a cubic volume V = R 3 . This is 
a simple approach as compared to more sophisticated ap¬ 
proaches needed to compute the proper density of states 
as discussed for various systems in Ref. mm- In principle 
there should be an infinite sum over discrete momentum 
values. However for simplification, we consider integra¬ 


tion over continuous values of momentum. Alternatively, 
surface and curvature effects should have been included 
in the continuous momentum variables which are beyond 
the scope of the present work. The essential physical ef¬ 
fects are expected to remain similar. Most likely we are 
going to make an underestimation of finite size effects in 
the quantities we evaluate. We have also not considered 
any modification of the model parameters with the expec¬ 
tation that much of those effects would show up in the 
mean values of the quantum fields. This is exactly in the 
line of neglecting the temperature and chemical potential 
dependence of model parameters in most of the related lit¬ 
erature. Within this approximation, the thermodynamic 
potential is now given by Eliza, 
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Here $ and are the Polyakov loop field and its conjugate 
respectively. The condensate fields of different flavors (/ = 
u, d, s) are given by a f = (ipftpf) = f x -fj d P > where 

E£ = yjp 2 + Mj is the corresponding quasiparticle energy 
with quasiparticle mass Mf = to/ — gs&f + gD&g&h- The 
NJL parameters are the 4-quark coupling gs, the 6-quark 
coupling gu and the ultraviolet cut-off A. The Polyakov 
loop potential is given by, 




= + <L 3 ) + 

—k ln[l - 6$i> + 4(4> 3 + <l 3 ) - 3($l>) 2 ] (2) 


where, b 2 {T) = a 0 + ai(-^) +a 2 (-p-) + a 3 (-+■) • Here a 0 , 
ai, 02 , To, ^3 and 64 are constants which may be obtained 
by fitting the temperature dependence of the Polyakov 
loop and pressure in pure gauge theory in the lattice 
framework (49]. The thermodynamic potential is first 
extremized with respect to the a and $ fields to estimate 
the mean field values at desired temperature and chemical 
potential. The values of the mean fields are plugged back 
into H(T, p) to obtain the mean thermodynamic potential. 

The model parameters used here are [32], ao = 6.75, 
oi = -1.95, o 2 = 2.625, o 3 = -7.44, b 3 = 0.75, 64 = 7.5, 
To = 190 MeV, k = 0.13, m u = uid = 5.5 MeV, m s = 
134.76 MeV, A = 631 MeV, g s A 2 = 3.67 and g D A 5 = 9.33. 
The cross-over temperatures for the different system sizes 
studied here is given in table [Q 

Results. — The globally conserved charges in the 
2+1 flavor matter are the quark number q, the electric 
charge Q and the strangeness S. The fluctuations of these 
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Fig. 1: (color online) Variation of second and fourth order susceptibilities as a function of T/T c for quark (left column), electric 
charge (middle column) and strangeness (right column). The volume dependence of c\ has been plotted in the inset. 


R(fm) 

2 

3 

4 

OO 

T c (MeV) 

160 

174 

178 

181 


Table 1: Cross-over temperatures corresponding to various sys¬ 
tem sizes as extracted in 2+1-flavored PNJL model at vanish¬ 
ing chemical potential. 


charges are related to the respective susceptibilities via 
the fluctuation-dissipation theorem, and are obtained as 
the moments of different orders from the corresponding 
chemical potential dependence of the free energy. Close 
to zero chemical potential these are given by, 


„x 


_ 1 d n (—VL(T, fi x )/T 4 ) , 

( ) n! d{^) n Imx- ° ( ) 


where, X stands for either of q , Q or S. These susceptibil¬ 
ities have been computed in first principle QCD calcula¬ 
tions on the lattice 17511531 as well as in calculations with 
hard thermal loops fHTT](TTT1 for large system sizes. Suit¬ 
able estimates of these fluctuations are also given by var¬ 
ious QCD inspired models (see e.g. [5TlfH51[fITll69llfi21f§3] ). 
In the context of finite volumes a study of fluctuations 
was done by some of us in Ref. [331 . Here we present the 
corresponding results in a 2+1 flavor PNJL model. 

To evaluate the fluctuations in the PNJL model, we first 
compute the thermodynamic potential D at any particu¬ 
lar temperature and with chemical potentials in a range of 
-300 MeV< fix <300 MeV with an interval of O.lMeV. 
The scaled thermodynamic potential ^, is then expanded 
in a Taylor series around fix/T= 0. The fluctuations of 
various orders may be extracted directly from the coeffi¬ 
cients of the series expansion. We have in fact fitted the 


scaled thermodynamic potential with a truncated Taylor 
series using the GNUPLOT software m- To obtain the 
best fit parameters up to the 4 th order, we have checked 
the least squares by varying the order of polynomials from 
6 to 10 as well as varying the range of fix/T around zero. 
The final least squares for the lowest system size were less 
than 10 -9 and much smaller (~ 10~ n ) for R=3fm, 4fm 
and oo. The whole procedure was repeated for the differ¬ 
ent temperatures in the range of 0.5 < < 2.5. 

In fig. □ we present the 2 nd and 4 th order susceptibil¬ 
ities for all the conserved charges as a function of t£-. 
From the inset we observe the saturation of Cj with in¬ 
creasing system size around T c . As in the case of 2 flavor 
system [33], the qualitative behavior of the susceptibil¬ 
ities for small system size is similar to that for infinite 
volumes, though quantitatively there are significant finite 
size effects. These effects are clearly visible for the 2 nd 
order susceptibilities in the range of 0.8T C < T < 2T c , 
and for the 4 th order susceptibilities in the range of 
0.8T C < T < 1.2 Tc. Therefore if the experimentally pro¬ 
duced fireballs thermalize and cool down to a temperature 
close to T c , one should consider freeze-out volume as a fit¬ 
ting parameter along with the temperature and various 
chemical potentials to fit the freeze-out multiplicity ratios 
(see e.g. m and references therein). Within the ambit 
of the present model, we find that the finite volume effects 
would be important if the system freezes out with a size 
less than R ~ 5/m. 

To find the relevant length scales in the different tem¬ 
perature regimes we have also checked the scaling of the 
susceptibilities with the corresponding mass scales in the 
model. For T <T C the relevant length scale is the temper- 
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Fig. 2: (color online) Scaling of susceptibilities with m„ below 
T c and with M q above T c as a function of T /T c . 


Fig. 4: (color online) Ratios of various correlations with fluc¬ 
tuations as a function of temperature. 


ature dependent pion mass ( m^R) and for T > T c the cor¬ 
responding scale is the temperature dependent constituent 
quark mass ( M q R , where M q is the average constituent 
mass for the three flavors). For demonstration we have 
shown these scaling behavior in fig. [2] The relevant masses 
were reported by some of us in |32j . We find that for very 
low temperatures as well as very high temperatures the 
mass scaling is satisfied. But in a significant temperature 
range close to T c the scaling is violated. 




T(MeV) T(MeV) 



T(MeV) 

Fig. 3: (color online) Ratio of fourth order to second order sus¬ 
ceptibilities for quark, electric charge and strangeness chemical 
potentials as a function of T/T c . 


It should be noted here that for comparing theoreti¬ 
cal computations like in the thermal model |100j . lattice 
QCD [101] . etc. with experimental data, one usually con¬ 
siders ratios of thermodynamic observables. This helps in 
situations where the volume is not directly measured and 
the free energy is assumed to be simply proportional to 
the volume of the system and therefore scales out in the 
ratios. In fact it has been argued (1011 that in such a situ¬ 
ation the susceptibility ratios are expected to be not only 
spectrum independent but dependent only on the ratio of 
q /T. Therefore Lattice QCD results corroborates hadron 
resonance gas results very well for temperatures below T c . 
With this expectation, attempts to match experimentally 
obtained fluctuation ratios with theoretical predictions are 
also undertaken [ 1021 . On the other hand some of us have 
reported that a finite size system of hadron gas may show 
some departure from volume scaling essentially in the elec¬ 
tric charge sector [34]. Therefore it is important to check 
the extent of volume scaling in the PNJL model. In this 
regard what we find is that the susceptibilities themselves 
being derivatives of the free energy density are still volume 
dependent. Therefore the free energy is not necessarily 
proportional to the volume of the system. So a corre¬ 
lated measurement of the particle multiplicities and the 
various susceptibilities may provide suitable estimates of 
the finite sizes of the produced strongly interacting matter 
with R < 5/m. Recently, with the assumption of volume 
scaling, an estimate of system size has been extracted from 
the fluctuation and correlations measured by the ALICE 
collaboration at y/s = 2.76 TeV confronting with the cor¬ 
responding computations in Lattice QCD 103 . 

On the other hand if the multiplicity ratios and ratios of 
fluctuations are to be employed to make contact between 
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Fig. 5: (color online) Ratios of net quark number density to 
potential at different temperatures. 

theory and experiment then it is important to find out if 
the volume scaling violation occurs for such ratios. Other¬ 
wise no unique value of the system volume may be found. 
In fig. J3|) we have plotted the ratios of fourth to sec¬ 
ond order fluctuations for different conserved charges as a 
function of temperature. It may be observed that the scal¬ 
ing of the ratios with the system volume holds good away 
from the crossover region. For all the conserved charges 
there is a significant violation of the volume scaling in a 
finite window of temperature. For the quark and electric 
charge sector the magnitude of violation is high compared 
to the strangeness sector, which however shows the scaling 
violation for the largest temperature range. 

Let us now carry out a study on ratios of leading or¬ 
der correlations to fluctuations among different conserved 
charges. The correlations are obtained as, 

X , F = 1 

<J iy! [ 1 

where, A" and Y stands for q , Q and S with X ^ Y. 
We compute the correlation coefficients using numerical 
differentiation. We obtain the thermodynamic potential 

at any given temperature within the range -50MeV< 
HXiHy <50MeV in an interval of 0.1 MeV. For obtaining 
the derivative g(f x /T)d(jj,J/T) ’ we com P u t e it for various 
values of A(/x x /T) = A /T) and consider the value 
when the derivative reaches a saturation. We repeat the 
whole procedure for the entire temperature range so as 
to obtain first order correlations of conserved charges as 
functions of temperature around zero chemical potentials. 

In fig. Q] we have plotted the ratios of correlations to 
fluctuations for different conserved charges. The volume 


the quark number fluctuations as functions of quark chemical 


scaling violation in this case is quite small for all the sec¬ 
tors except for the smallest system size. We note that 
here the correlations and fluctuations are derivatives of 
same order of the free energy. 

The fluctuations of the conserved charges can also be 
extracted for non-zero chemical potentials. Here we have 
considered net quark density c\ and quark number fluc¬ 
tuations C 2 for simplicity. These quantities are also the 
simplest to analyze from the experimental data at least 
in terms of net proton number and its fluctuations 1102 ] . 
We compute the thermodynamic potential as a function 
of fj, q in a range of 0 < fi q < 400MeV with an interval of 
O.lMeV for a set of temperatures. Then using the method 
of numerical derivative we obtain the observables. 

The quark number density itself was found to have a 
volume dependence. The ratio c\/c\ is displayed in fig. [5] 
for a few representative values of temperature. Since the 
quark number is very small for low values of /j, q we present 
our results in the range 100MeV< p , q <400MeV. For high 
temperatures the ratio is monotonically increasing. On 
the other hand for T < 200 MeV we observe a dip in a 
certain window of n q . Incidentally these windows of ^, q 
correspond to those close to the phase boundary [32] . The 
dips occur due to the shooting up of the quark number 
fluctuations in these regions. For T < 100 MeV there is 
even a discontinuity indicating the existence of the first 
order line in the phase diagram for the corresponding sys¬ 
tem size. We observe that for high temperatures the vol¬ 
ume scaling is maintained even for large chemical poten¬ 
tials. However as we go down below T ~ 200 MeV the 
volume scaling is violated in the transition region of the 
corresponding temperature. Thus we may infer that the 
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violation of volume scaling if any, always occur close to 
the transition region, where large correlation lengths come 
into play and lead to the separate finite size behavior of 
the derivatives of the free energy. 

Conclusion. — To conclude, we have studied the vol¬ 
ume dependence of the free energy density of 2+1 flavor 
strongly interacting matter in terms of various second and 
fourth order fluctuations of the conserved charges, as well 
as the quark number density at finite temperature and 
chemical potentials. We observed volume scaling viola¬ 
tions at two levels. Firstly the free energy density is itself 
volume dependent for R < 5/m within a certain range of 
T and n q , as given by the behavior of the susceptibilities. 
Secondly the ratio of these derivatives themselves show vi¬ 
olation of volume scaling in a small window of T and 
all along the transition region in the T — phase dia¬ 
gram. The addition of one more flavor in our treatment 
did not change any qualitative features of the quark num¬ 
ber fluctuations compared to the the 2 flavor system pT51 . 
though the quantitative results are significantly different. 
The 2nd order strangeness susceptibility has a substan¬ 
tial contribution to the quark number fluctuation. On the 
other hand the 4th order strangeness susceptibility has 
a sub-dominant contibution to the corresponding quark 
number susceptibility. As a result the volume scaling vi¬ 
olation in the ratio of 4th order to 2nd order quark num¬ 
ber fluctuations is somewhat smaller in 2+1 flavor system 
compared to the 2 flavor system. However a cleaner signal 
of volume scaling violation is now observed in the ratios of 
strangeness susceptibilities themselves, which is sustained 
in a much wider range of temperatures. We have also been 
able to extract various charge correlations. The ratios of 
these correlations with the 2nd order fluctuations violate 
volume scaling only for very small system sizes. We have 
further extended our studies to non-zero chemical poten¬ 
tials. An interesting feature observed is that the first or¬ 
der phase transition nature is itself manifestly dependent 
on the system size, as can be seen from the quark num¬ 
ber to quark susceptibility ratios. For very small systems 
the first order phase boundary is completely washed out. 
Further investigation in these directions are required with 
the computation of higher order moments that we wish to 
carry out in future. 
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